Nonlocal field correlators on the lattice in HP cr-model 
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Abstract 

Connected two-point field strengtli correlators have been measured on the lattice in quatcr- 
nionic projective cr-model of pure SU(2) Yang-Mills theory. The correlation lengths, ex- 
tracted from the exponential fit for these correlators, are found to be AJ"^ = 1.40(3) GeV and 
A^^ = 1.51(3) GeV in good agreement with other existing calculations. The dependence of 
bilocal functions on the connector shape was studied. 
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The studies of nonperturbative (NP) aspects of QCD are well known to be of prime 
importance. This is the most interesting and rich field, albeit the most complicated one. It 
is commonly believed that the origin of complex NP phenomena such as confinement, chiral 
symmetry breaking etc is highly nontrivial structure of QCD vacuum. The important step in 
the NP vacuum fields investigations was undertaken in the seminal paper [1] where the gluon 
condensate - NP average of gluon fields over vacuum state - was introduced. This object, 
still playing essential role in studies of QCD vacuum structure is defined as 

f)/ 

G2 = ^{FiMKM). (1) 

where F'^^^'s, are the field strength tensors. In the QCD sum rules approach this is universal 
quantity (together with chiral condensate (V'V') and a few higher NP averages), describing 
QCD vacuum. The sum rule method and idea about NP condensates have proven to be very 
fruitful. Using this technique masses of resonances, decay constants and other quantities 
of phenomenological interest were computed in good agreement with experiment. However, 
there are questions, which are difficult to answer in the sum rule method. In general, they 
arise when one is interested in such objects as potentials, excited states etc, i.e. for physical 
situations when one of the space-time scales characterizing the problem in question is large. 
It was found (in [2' on the lattice and later by other authors), that the word "large" here 
has rather precise meaning: there is an important parameter, characterizing NP dynamics 
of vacuum fields - correlation length, which defines the spacial decay of connected gauge- 
invariant bifocal correlators of field tensors. It is natural to expect that original sum rule 
method, operating with local quantities, is not applicable to the situations where nonlocal 
properties of vacuum fields are essential. 

The method of field correlators (MFC) |3l 01 15] (see j6] for review) can be considered as 
a development of sum rules method explicitly taking into account long distance effects. The 
dynamical input is a set of NP field strength correlators 

^i,^v^,...,i,^v^ =1^ {,gF^,^y^{xx)^{xx,x2)gF^^y^{x2) ■ . .gF^^y^{xn)^{xn,x{)) , (2) 
/. \ 

where ^>(x,y) = Pexp J A^dz^ with integration along some path, connecting the points 

\y J 

X and y, are the phase factors, introduced for the sake of gauge invariance. As a matter of 
principle these objects can be understood as infinite power series of NP condensates and in 
this sense their theoretical status is absolutely equivalent to ([T]). But it is more convenient 
as is seen below, to express ([2]) in terms of correlation lengths. Gauge-invariant observables 
for QCD processes of phenomenological interest can be expressed through ^ via cluster 
expansion. MFC includes not only perturbative QCD at small distances, but also large 
distance nonperturbative effects (see |6j and references therein). 

It is remarkable, that one can describe most NP QCD phenomena with high accuracy by 

(2) 

means of the lowest two point correlator while higher cumulants can be considered 

as small corrections [6J. In particular, the direct consequence of this fact is the so called 
Casimir scaling law for static potential between two quarks [71 [8] , clearly seen on the lattice 
[9]. The corresponding formalism is called the Gaussian dominance approximation. Here we 
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have single fundamental input - bilocal correlator, which has the following parametrization 



d 



where L'(x^) and Di{x'^) are some functions of distance between two points and integration 
in <I>(x,0) goes along the straight lines. 

At short distances the expression ^ can be expanded in powers of x"^, the corresponding 
coefficients happen to be proportional the condensates of higher powers (see [TOl [TT| [T2| [T3j [Ti] 
for details). The leading nonabelian term is given by (notice that our normalization is different 
from that of the cited papers): 

dD{x) _ / rabc,jpa j:,b 



dx"^ 



rabc, pa pb pc \ 
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(thus at small distances the correlator is Gaussian as a function of x). On the other hand 
the NP long distance contribution can be parameterized by exponentials with the correlation 
lengthes A and Ai: 

D{x) = Aexp(-|x|/A), Di{x) = exp(-|2;|/Ai). (5) 

Our lattice results demonstrate that exponential behavior is a good approximation starting 
just from the distance of order of the correlation length A, Ai (see details below), while for 
smaller distances the regime ^ is expected (but not explicitly seen by us due to insufficient 
lattice precision). 

Very important source of information on field correlators are lattice measurements [21 
[T5] . There is an immediate problem, however. Measuring ([s]) on the lattice one would get 
~ piece at small distances as predicted by perturbation theory together with some 

NP contributions like ^ and extracting the NP part is a nontrivial task. This problem 
can be solved by using cooling technique (as in [2]), which allows to eliminate short range 
fluctuations, or by smearing procedure as in [15]. Another possibility is to expand nonlocal 
correlators in powers of local condensates, but needless to say that practically it is not possible 
to study all the series of local quantities. So, it is desirable to find some alternative approach, 
which allows to single out the NP signal and to compare with the other lattice methods. 

For this purpose we can use an opportunity to consider modified field configurations, 
close in some sense to the initial configurations. A nice example of such modification is re- 
placement of the initial SU(2) gauge theory by configurations of scalar fields representing 
nonlinear a-model with target space being quaternionic projective space HP^. This was re- 
alized on the lattice in \16\ [T7] , where both local objects like gluon condensate and nonlocal 
quantities were measured in good agreement with known results. So, the string tension calcu- 
lated in terms of HP^ projected fields turned out to be very close to full SU(2) string tension. 
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\^a^^^ /a^^^'^^Y^'^ = 1.04(3) in the continuum limit. It was established, that this projection 
captures only NP content of gauge background, cutting contribution from perturbation the- 
ory. So, it is very natural to look at correlators of gluon fields in this approach and this is 
done for the first time in the present paper. 

Due to lack of space we can not discuss all details of constructing HP-projected fields (see 
|16] for this purpose), but to be self-contained we give here the main steps. The simplest 
explicit parametrization of HP^ is provided by normalized quaternionic vectors 

k> = Qi G H, {q\q) = q^qi = 1 e H, (6) 

where H is the field of real quaternions. The states \q) describe 7-dimensional sphere, while 
the HP^ space is the set of equivivalence classes of \q) with respect to the right multiplication 
by unit quaternions (elements of SU(2) group) 

\q)r^\q)v, \vf = l, v€H. (7) 

Configuration \qx), assigned to the lattice as the best possible HP^ fields approximation 
provides the minimum of functional 



F{A,q) = j{A, + {q\d,\q)Y (8) 
for given SU(2) field A^. The working tool for computations is link variable 

TT _ {Q.x\Q.X + fi) /„X 

\{(lx\qx+^Ji)\ 

What is actually measured is a loop made of two plaquettes (smallest rectangular contours 
made of four link variables) P^xuix) = Ux,f_iUx+fj.,uU^^^^^U^,^ connected by the phase factors 
(product of links Ui between two points) along straight lines, with different orientation of the 
plaquettes: 



a 




^ (tt{p^,,,{x) -i)[[[u,] {P,,uM -i)[YYut\) , (10) 




where a is the lattice spacing. 

There are two possibilities to obtain nonzero answer in (|3]) - when parallel planes {jiivi) 
and (/i2i^2) are perpendicular to the vector x {D±{x)) and when planes (/^iz^i) and {(^21^2) are 
parallel and x lies in this plane (D^^{x)): 

D±{x) = D{x)+Di{x), (11) 
Z?l,(x) = L»(x) + Di(x)+x2^^iM. (12) 

One attractive feature this approach is its relatively low computational cost. Lattice 
simulations have been performed on PC 1.66GHz on 100 configurations (taken from |18j ) 
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at /3 = 2.6 corresponding to the lattice spacing a = 0.06 fm. The total volume was 40^. 
Other configurations with larger lattice spacing were also used for scaling properties study. 
Parametrization ([5]), which works very well starting from two lattice spacing (remember that 
the HP^ projected fields do not contain perturbation theory contributions at small distances) 



gives for correlator Di{x) from (11) (see Fig.l) rather stable result A^^ = 1.40(3) GeV (or 
Ai = 0.14(1) Fm) in good agreement with [19]. The data for the correlator D(x) (Fig. 2) 
allow us to extract A~^ = 1.51(3) GeV (or A = 0.13(1) Fm). The preexponential factors 
A « 0, 11 GeV^, Ai 0,06 GeV^. We studied the scaling properties of correlation length and 
found that the dependence on the lattice spacing is almost inessential for this quantity (see 
Fig. 3). On the other hand, the quantity A + Ai, which is proportional to the condensate G2, 
scales as a/a^ + /3 (see Fig.4) with a ^ (95MeV)2 and G2 = {6Nj7r'^)P = 0,062(6) GeV^ in 
complete agreement with results of |17] . The quoted uncertainties are the standard statistical 
errors. For reader's convenience we give the values of gluonic correlation length A, obtained 
in other papers, in Table 1. 

Table 1. Correlation length A. 



U3\SU{2) [2,SC/(3) |l5]5f/(3) W\SU{3) 

A, fm 0.13 0.22 0.12 0.10 



It is interesting to check how the correlators depend on the shape of curve connecting 
the points and x. It turned out, that correlator D±{x) holds the form close to exponential 
decay, if we replace straight connector by the U-like curve (see Fig. 5). Fig. 6 shows how that 
exponential index A depends on the depth / of the curve flexure. Unlike to this case of sizeable 
changing of correlation length for the function D±(x) (by ~ 16%) correlator -D||(x) does not 
preserve exponential form under the changing of the shape. We have also checked how shape 
of the connector influences the integral quantities. To this end, we consider the sum of 
two-point "U-like" and "L-like" correlators, taken over rectangle R x T in the (xt)-plane 

1,1' 

The typical contributions to the sum are shown on Fig. 7. The value wu^l scales as the 
rectangle square wu^l «; ajj^iRT for sufficiently large rectangles. The perimeter terms, also 
presented in wjj^l-, may be taken into account as R- and T-dependence of au,L- The depen- 
dence of correlators on the connector shape shows as difference between au and (Tl ■ The sum 



(13) as function of rectangle square is shown on Fig. 8 at fixed R = 7. As turned out the 
difference is small: 

^ 5%, (14) 

so the integral quantities are almost insensitive to changing of the connector form. Another 
interesting question is the relation of the coefficient ct[/,l to the string tension a^^^ (one 
would expect approximate equality of these quantities in Gaussian stochastic model). It will 
be discussed elsewhere 1211. 
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Finally, it is natural to expect that correlators of pure electric and pure magnetic fields 
8{x) = Ag]4(x), B{x) = leumei.k^fl^ki^) (15) 



should coincide at x = because of 0(4) rotational symmetry of the problem. As we can see 
from Fig. 9, it is indeed so - two lines, corresponding to logarithms of electric and magnetic 



correlators (15) come together at small x (the error bars are within the symbols). 

In conclusion, we have measured the two-points correlators of vacuum fields on the lat- 
tice in HP^ (T-model. This approach has allowed us to get rid of the perturbation theory 
contribution and to use the long-distance nonperturbative expression in the whole region of 
distances, beginning from two lattice spacing. Our measurements provide solid confirma- 
tion of exponential behavior for long distances of functions D{x) and Z)i(x), parameterizing 
two-point correlators. Small distance behavior is the question of separate interest, but due 
to insufficient precision of our lattice calculations we can not catch expected Gaussian x^- 
dependence and we do not concentrate on this problem, being interested in long-distance 
characteristics of two-point correlators. So, the corresponding correlation lengthes are found 
to be A^^ = 1.40(3) GeV and = 1.51(3) GeV in good agreement with other known cal- 
culations using completely different lattice technique. It is worth noticing that parameters 
^ = {AX^Y/'^ 0.14, ^1 = (^iAf)^/2 QJ2 turn out to be rather small, confirming the MFC 
picture of NP vacuum fields being weak in units of the inverse correlation length (see related 
discussions in [22]). We also studied the dependence of bilocal functions on the connector 
shape and found some decrease of the correlation length with the increasing of the connector 
length, corresponding to the function D±{x). In the case of the function D\\{x) this depen- 
dence can be established via the integral characteristic as the small difference (about 5%) 



between slopes of wu{S) and wl{S) (13) as functions of rectangle square. 

The authors are grateful to F. V. Gubarev, M. I. Polikarpov and Yu. A. Simonov for helpful 
discussions. The research was supported by the grant for support of scientific schools NS- 
4961.2008.2. 

References 

[1] M. A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B147, 385; 448 (1979). 
[2] A. Di Giacomo, H. Panagopoulos, Phys. Lett. B285, 133 (1992). 
[3] H.G. Dosch and Yu.A. Simonov, Phys. Lett. B205, 339 (1988). 
[4] H.G. Dosch, Phys. Lett. B190, 177 (1987). 
[5] Yu.A. Simonov, Nucl. Phys. B307, 512 (1988). 

[6] A. Di Giacomo, H.G. Dosch, V.I. Shevchenko, Yu.A. Simonov, Phys. Rept. 372, 319 
(2002), |hep-ph70007223[ 



[7] Yu.A. Simonov, JETP Lett. 71, 127 (2000), hep-ph/000124^ 



[8] V.I. Shevchenko, Yu.A. Simonov, Phys. Rev. Lett. 85, 1811 (2000), hep-ph/0001299 



6 



G.S. Bali, Nucl. Phys. Proc. Suppl. 83, 422 (2000), |h ep4at/9 90802T| 
S.N. Nikolaev and A. V. Radyushkin, Nucl. Phys. B213, 285 (1983). 
S.V. Mikhailov, Phys. Atom. Nucl. 56, 650 (1993). 
A.G. Grozin, Int. J. Mod. Phys. AlO, 3497 (1995), [hep^ph/9412238[ 



A.E. Dorokhov, S.V. Esaibegian and S.V. Mikhailov, Phys. Rev.D56, 4062 (1997), 
ppjph/97024lf[ 

V.I. Shevchenko, Yu.A. Simonov, Phys.Atom.Nucl.60, 1201 (1997), [hep^th/9701026 



G.S. Bah, N. Brambilla, A. Vairo, Phys. Lett. B421, 265 (1998), pIeFlat/9709079[ 
F.V. Gubarev, S.M. Morozov, Phys. Rev. D72, 076008 (2005), |hep-lat/0509011[ 
P.Yu. Boyko, F.V. Gubarev, S.M. Morozov, Phys. Rev. D73, 014512 (2006), ,hep^ 



lat/0511050| 



F. V. Gubarev, private communication. 

M. Campostrini, A. Di Giacomo, G. Mussardo, Z. Phys. C25, 173 (1984). 

A.M. Badalian, A.V. Nefediev, Yu.A. Simonov, JETP Lett. 88, 558 (2009), 



arXiv:0807.3320 [hep-ph]. 



V.D. Orlovsky, V.I. Shevchenko, in preparation. 

Yu.A. Simonov, V.I. Shevchenko, |arXiv:0902.1405' [hep-ph]. 



7 




Figure 1: Function logZ^i(a;) (conventional units) from the measurements of two point correlators 
(9) 




Figure 2: Function log D{x) (conventional units) from the measurements of two point correlators 
(9) 
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Figure 3: Correlation length Ai as a function of lattice spacing 




Figure 4: Condensate -D(O) + -Di(O) as a function of lattice spacing 
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Figure 5: Correlator on the lattice for U-like connector 
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Figure 6: Exponential index as a function of the connector shape 
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Figure 7: Typical contribution to the sum (12) for (a) "U"-like and (b) "L"-like correlators 
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Figure 8: The sum (12) as function of rectangle square S = {R = 7) x T ior "U"- and "L"-like 
correlators 
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Figure 9: Correlators of electric and magnetic fields (14) 
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